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ABSTRACT 

Individual pulses from pulsars have intensity-phase profiles that differ widely from pulse to 
pulse, from the average profile, and from phase to phase within a pulse. Widely accepted expla- 
nations for pulsar radio emission and its time variability do not exist. Here, by analysing data 
near the peak of the Vela pulsar's average profile, we show that Vela's variability corresponds to 
lognormal field statistics, consistent with the prediction of stochastic growth theory (SGT) for 
a purely linear system close to marginal stability. Vela's variability is therefore a direct mani- 
festation of an SGT state and the field statistics constrain the emission mechanism to be linear 
(either direct or indirect), ruling out nonlinear mechanisms like wave collapse. Field statistics are 
thus a powerful, potentially widely applicable tool for understanding variability and constraining 
mechanisms and source characteristics of coherent astrophysical and space emissions. 

Subject headings: methods: statistical — plasmas — pulsars: general — pulsars: individual (Vela) — 
radiation processes: non-thermal — waves 



1. Introduction 

Pulsars are highly magnetized neutron stars 
whose rotation causes highly nonthermal beams 
of radiation to be swept across the Earth (Manch- 
ester & Taylor 1977), similar to the periodic view- 
ing of a lighthouse beam. Most likely the radio 
emission is produced over the star's (magnetic) 
polar caps (Manchester & Taylor 1977; Melrose 
1996; Asseo 1996). Since their discovery in 1967 it 
has been recognized that only suitably long time 
averaging leads to a stable intensity profile versus 
pulsar phase. While this average profile is unique, 
individual pulses vary widely in intensity, often 
by a factor of 5 or more, from one phase to an- 
other in a given pulse and from one pulse to the 
next at a given phase, as in Fig. 1. This vari- 
ability (Manchester & Taylor 1977; Hankins 1996) 
includes phenomena known as drifting sub-pulses 
(Drake & Craft 1968), microstructure (Craft et al. 
1968), giant pulses (Cognard et al. 1996) and giant 
micropulses (Johnston et al. 2001). Sub-pulses are 
features which drift in time across the pulse win- 
dow whereas microstructure are concentrated fea- 
tures superposed on a subpulse that are sometimes 



quasiperiodic. Giant pulses and micropulses are 
very rare pulses with fiuxes ~ 10 times the average 
flux (Cognard et al. 1996; Johnston et al. 2001). 
No accepted explanation exists for these forms of 
variability or, indeed, for the mechanism(s) pro- 
ducing pulsar radio emission (Manchester & Tay- 
lor 1977; Mehose 1996; Asseo 1996; Hankins 1996). 

The high brightness temperatures of pulsars re- 
quire coherent emission processes such as plasma 
microinstabilities or nonlinear processes. Linear 
mechanisms include (Melrose 1996; Asseo 1996; 
Luo & Melrose 1995):(1) linear acceleration and 
maser curvature emission, in which electrons radi- 
ate coherently while accelerating in an oscillating 
large-scale field or on curved magnetic field lines, 
respectively, and (2) relativistic plasma emission, 
in which a streaming instability either directly 
generates escaping radiation near harmonics of the 
electron plasma frequency fpe or else drives local- 
ized (non-escaping) waves near fpe that are con- 
verted into escaping harmonic radiation by linear 
mode conversion or nonlinear processes. Nonlin- 
ear mechanisms can produce radiation by wave co- 
alescence and scattering processes or as intense lo- 
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calized wavepackcts, perhaps driven near /pe by a 
streaming instability, undergo modulational insta- 
bilities and strong turbulence wave collapse (Asseo 
et al. 1990; Asseo 1996; Weathcrall 1998). Since 
existing analyses suggest that many mechanisms 
are viable, in part due to large uncertainties in 
the plasma properties and location of the emit- 
ting regions (e.g., above the polar cap or near the 
light cylinder), new approaches are necessary. 

Analyses of intensity scintillations and angu- 
lar broadening, corresponding primarily to Fourier 
analyses of data, are standard for astrophysical 
and solar system radiation sources (Rickett 1990). 
In contrast, distributions of electric field strengths 
or intensities were rarely analyzed until recently, 
perhaps because their strong theoretical motiva- 
tions and benefits were not clear before the ad- 
vent of stochastic growth theory (SGT) (Robin- 
son 1992; Robinson et al. 1993; Cairns and Robin- 
son 1999; Cairns et al. 2000; Robinson and Cairns 
2001) and other theories like self-organized criti- 
cality (SOC) (Bak et al. 1987). However, recent 
analyses of 7 different solar system wave phenom- 
ena show that all have well-defined field distribu- 
tions that agree very well with the predictions of 
SGT (Robinson et al. 1993; Cairns and Robin- 
son 1999; Cairns et al. 2000; Robinson and Cairns 
2001), resolving longstanding theoretical problems 
pertaining to the burstyness, widely varying fields, 
and persistence of the waves. Similarly, the giant 
pulses of some pulsars have power-law fiux distri- 
butions (Cognard et al. 1996), sometimes inter- 
preted qualitatively in terms of SOC (Young & 
Kenny 1996). With the advent of rapid time res- 
olution, coherently de-dispersed data for Vela and 
other pulsars (Johnston et al. 2001) the time is 
ripe for analyzing pulsar variability and its statis- 
tics in terms of SGT, SOC and other theories. 

This paper directly addresses pulsar variability 
and emission mechanisms by analysing the radi- 
ation's statistics near the peak of the Vela pul- 
sar's average profile and interpreting the results 
in terms of the theoretical predictions and formal- 
ism of SGT. After summarizing the predictions of 
SGT and other theories for wave growth (Section 
2), Section 3 shows that Vela's intrinsic variability 
near the peak of the average pulse profile corre- 
sponds to lognormal statistics in the electric field 
(or intensity), not Gaussian or power-law statistics 
in the intensity. Pulse variability is thus a direct 



manifestation of an SGT state. The consistency 
with the SGT prediction then strongly constrains 
the emission mechanism and source plasma (Sec- 
tion 4) , with nonlinear emission mechanisms being 
non-viable in the phase range analysed. Prelimi- 
nary results at other phases of Vela's pulse profile 
and for other pulsars are then briefly discussed 
(Section 5). This analysis provides a first demon- 
stration that radiation statistics for astrophysi- 
cal sources are a powerful and potentially widely 
applicable tool for strongly constraining emission 
mechanisms and source plasmas. 

2. Theories for Wave Statistics 

Wave-particle interactions are expected to drive 
natural plasmas towards marginal stability, where 
wave emission and damping (as well as total en- 
ergy inflow and outflow) are balanced. SGT treats 
systems in which an unstable particle distribution 
interacts self-consistently with its driven waves in 
an inhomogeneous plasma background and evolves 
to a state in which (i) the particle distribution 
is close to time- and volume-averaged marginal 
stability but with stochastic fluctuations that (ii) 
cause the wave gain G to be a stochastic variable 
(Robinson 1992; Robinson et al. 1993; Cairns and 
Robinson 1999; Cairns et al. 2000; Robinson and 
Cairns 2001). Here G{t) is related to the wave 
growth rate T{t) by G{t) = J^^ dt'T{t') and to 
the time- varying wave electric field E{t) and a ref- 
erence field Ea by -E^(t) = Eq exp G{t). Rewriting 
this time integral as a summation over fluctuations 
AGi = TjJ^ti then, provided only that sufficiently 
many fluctuations in AGi occur in some charac- 
teristic time, the Central Limit Theorem requires 
that G{t) is a Gaussian random variable irrespec- 
tive of the detailed distribution of AG^. Hypothe- 
ses (i) and (ii) thus have simple and natural phys- 
ical justifications. The hypothesized random walk 
in G oc ln£^ then implies that the waves should 
be bursty and widely varying in amplitude, while 
the closeness to marginal stability implies that 
the waves and driving distribution should persist 
far from the latter's source. These characteristics 
are very attractive for pulsars, given the existence 
of intrinsic variability and the radiation's broad 
bandwidth (and so large radial extension of the 
source inferred therefrom), as well for many other 
astrophysical and space phenomena. 
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Due to G being a Gaussian random variable, 
pure SGT predicts that the probabihty distribu- 
tions of wave field and intensity are lognormal 
(Robinson ct al. 1993; Cairns and Robinson 1999; 
Robinson and Cairns 2001); i.e., 

P{\ogE) = (\/2^ a)-^ exp(-(logi; - nf I2a^) . 

(1) 

where log means to the base 10, /x and a are the 
average and standard deviation of logi?, respec- 
tively, and / d{[ogE)P{[ogE) = 1. Nonlinear 3- 
wave processes active at high E above a thresh- 
old Ec, which remove energy from the waves, re- 
duce the P{\.ogE) distribution below the predic- 
tion (1) near and above E^ with known analytic 
form (Robinson ct al. 1993; Robinson and Cairns 
2001). Processes like wave collapse and modu- 
lational instability cause a power-law tail with 
P{E) cx E~°' , with a ranging from 4 to 6, to 
develop above E^ (Robinson 1997; Robinson and 
Cairns 2001). Waves driven from thermal levels by 
an instability, and which retain memory of their 
thermal past, also develop a power-law tail, but 
usually with a smaller index (Cairns et al. 2000; 
Robinson and Cairns 2001). Finally, SOC should 
produce a power-law distribution with index close 
to —1 (Bak et al. 1987) and the usual model for 
wave growth in plasmas (uniform secular growth 
with constant growth rate) should produce a uni- 
form distribution at fields below Ec (Robinson et 
al. 1993; Cairns and Robinson 1999; Robinson and 
Cairns 2001). In contrast, scattering by density 
turbulence or radiation from multiple incoherently 
superposed sources is expected to produce Gaus- 
sian intensity distributions (Rickett 1990). Rig- 
orous testing of theories for wave growth is thus 
possible using the observed field statistics, as al- 
ready demonstrated in multiple space contexts ref- 
erenced above. 

3. Statistics of Vela's intrinsic variability 

The data set consists of 20085 contiguous pulses 
(30 minutes) of the Vela pulsar, measured at 
1413 MHz by the Parkcs radio telescope and pro- 
cessed using coherent dedispersion and other tech- 
niques by Johnston et al. (2001). There are 
2048 phase (time) bins per pulse period, each of 
44 /xs length (comparable to the scatter broaden- 
ing time). Vela's average intensity (over many 
pulses) is restored by adding Iq = 1250 mJy to 



each sample. Fig. 1 shows the average pulse pro- 
file for relevant phase bins in mJy, together with 3 
superposed pulses which illustrate the variability. 
Note that the noise level is very low compared with 
earlier analyses, allowing detailed investigation of 
the intrinsic field statistics. 

Analysis of data in the off-pulse phase bins 
lead to Gaussian statistics in the intensity I, 
as expected for instrumental and background 
noise. For instance, fitting the P{I) distribu- 
tion for phases 391-399 to a Gaussian, by using 
the Amoeba algorithm to minimize (Press et 
al. 1986), yields (/) = 1215 m.Jy (agreeing with 
Iq to within less than the 100 mJy bin width), 
a I = 1420 mJy, = 66 for Ndf = 46 de- 
grees of freedom, and a significance probability 
-P(X^) = 0.03. (Fitting is restricted to intensity 
bins with > 100 pulse samples). This fit has good 
statistical significance. 

Fig. 2 shows the P{I) distribution and its best 
Gaussian fit for phase bin 490, close to the peak 
in Vela's average profile. The fit clearly fails at 
both low and high /, entirely missing the long tail 
at large /, as confirmed by it having = 301 for 
Ndf = 53 and P{x^) < IQ-^e. The variability 
at this phase is thus not described by Gaussian 
intensity statistics. 

In contrast, defining the electric field E = 
/(mJy)i/2, Fig. 3 shows that the P{\ogE) dis- 
tribution for phase 490 is well fitted by the SGT 
prediction (1): for bins with > 100 pulse samples 
and E > 10^ imits (intensities above 10"* mJy, 
which is 6(7/ above the noise), /U = 2.3, a = 0.096, 
= 27 for Ndf = 19, and P(x^) = 0.12. 
The Kolmogorov-Smirnov test (Press et al. 1986) 
yields a significance probability of 47%. This fit 
is strongly statistically significant, clearly demon- 
strating that pulsar variability at this phase is 
lognormally distributed and quantitatively consis- 
tent with the theoretical form predicted by simple 
SGT. The fit matches the data well even outside 
the fitted range of fields (dotted line), although 
the effects of the noise background become increas- 
ingly evident at fields < 80 units. 

Results similar to Fig. 3 are found for phases 
485 - 540, for which the average pulsar intensity is 
well above the noise level, although the statistical 
significance varies. Rather than showing more re- 
sults for individual phases. Fig. 4 shows the P{X) 
distribution observed for phases 485 - 500 simul- 
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tancously, where X — (logi? — fi{4')) / a{(j)) is the 
field variable resulting from detrending variations 
in jj, and a with phase <p. Comparison with (1) 
shows that SGT predicts the P{X) distribution to 
be Gaussian with zero mean and unit standard de- 
viation (solid curve) (Cairns and Robinson 1999). 
The agreement is very good, with the Kolmogorov- 
Smirnov test yielding a significance probability of 
0.1%. 

4. Theoretical Implications for Vela 

For phases 485 - 540, where Vela's average pulse 
profile is well above the noise, the field distribu- 
tions do not have power-law tails or nonlinear cut- 
offs. Instead, the data have lognormal statistics 
and the variability is a direct manifestation of a 
simple SGT state, with no evidence for nonlin- 
ear processes, SOC, or uniform secular growth. 
This absence of a power-law tail or cutoff in the 
P(log£^) distributions for these phases rules out 
pulsar emission mechanisms based on nonlinear 
processes (Asseo et al. 1990; Asseo 1996; Weather- 
all 1998) such as wave collapse, modulational in- 
stability, and three-wave processes. Instead, the 
observed consistency with simple (linear) SGT 
means that only linear emission mechanisms arc 
viable, meaning that a plasma instability in a SGT 
state either directly generates the radiation or else 
generates non-escaping waves that are tranformed 
into escaping radiation by linear processes (e.g., 
mode conversion) alone. 

Prom the definitions of fj, and a and the in- 
tensity decreasing with distance R as R~^, it 
is easy to show that cr(i?) is constant and that 
/x(i?) = n{Ro) — log(i?/i?o); where ij,{Ro) is the 
value at the source's edge (R ~ Rq). Taking the 
values n = 2.0 and a = 0.1 to be representative 
of these phases, the distance R = 350 pc for Vela, 
and the value Rq = 30 m, yields /i(i?o) ~ 20. The 
value i?o = 30 m results from assuming that the 
overall source is annular, with radius equal to the 
neutron star radius « 10 km, and dividing by the 
2048 phase bins used for Vela. Accordingly, the 
ratio /txo/c « 200 in the source. The values /io and 
a will constrain future theoretical models for why 
SGT applies. 



5. Discussion and Conclusions 

The foregoing analyses are the first applications 
of SGT to propagating EM radiation and, simul- 
taneously, to extra-solar system sources. Their 
success implies that radiation statistics are an un- 
derappreciated and potentially very powerful tool 
in astrophysics (and space physics), and suggests 
that SGT may well be widely applicable to coher- 
ent astrophysical sources. As to whether the Vela 
results are representative of other pulsars, analyses 
arc ongoing. Our results to date for pulsar PSR 
1641-45, see also Johnston & Romani (2001), sug- 
gest that the variability near the peak of the av- 
erage profile also corresponds to lognormal statis- 
tics, thereby being consistent with SGT and the 
Vela results above. 

Of course, SGT is not likely applicable to all 
sources or indeed to all components of pulsar emis- 
sions. For instance, Jovian "S bursts" have a 
power-law flux distribution with index 2.0 ± 0.5 
(Queinncc & Zarka 2001) and the peak flux dis- 
tribution of solar microwave spikes can be fltted 
with an exponential or perhaps a lognormal form 
(Isliker & Bcnz 2001). Moreover, this richness 
in possible wave statistics also appears in phase 
bins away from the peak in Vela's average pro- 
file and for pulsars with giant pulses, where the 
observed P{logE) distributions are often approxi- 
mately power-law. Detailed interpretations will be 
described in detail elsewhere. For now, we men- 
tion only that indices « 4.5±1.0 are likely too high 
for SOC but are instead probably due to either 
driven thermal waves (Cairns et al. 2000) and/or 
due to strongly nonlinear processes like modu- 
lational instability and wave collapse (Robinson 
1997; Robinson and Cairns 2001). The latter idea 
complements earlier suggestions (Asseo et al. 1990; 
Weatherall 1998) and appears particularly attrac- 
tive for giant pTilses and giant micropulses. 

In conclusion, analysis of rapidly-sampled, co- 
herently de-dispersed data near the peak of the 
Vela pulsar's average intensity-phase profile show 
that the field statistics are lognormal and quan- 
titatively consistent with SGT's prediction for a 
purely linear system near marginal stability. The 
variability is thus a direct manifestation of an SGT 
state and only linear emission mechanisms (either 
direct or indirect) arc viable. Observations for 
other pulsars and at other phases for Vela yield 
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both similar and different results, hinting at a pos- 
sible richness of wave statistics and emission mech- 
anisms. Analysis of field statistics is thus a pow- 
erful tool for understanding source variability and 
constraining the emission mechanisms and source 
characteristics that may be widely useful for coher- 
ent astrophysical and solar system radio emissions, 
as already found for plasma waves in space. 
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450 500 550 

Phase Bin 

Fig. 1. Average intensity-phase profile of the 
Vela pulsar for the dataset (thick solid line) to- 
gether with three superposed individual pulses 
(other lines). 



ntensity (mJy) 

Fig. 2. — The observed distribution P{T) at phase 
490, formed by binning the data linearly in the in- 
tensity and normalizing by J dlP(I) = 1, is shown 
with open circles and ±VlV error bars. The solid 
line shows the best-fit Gaussian, calculated for the 
domain shown by the dashed line. 
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Fig. 3. The observed distribution P(logii') at 
phase 490 (open circles and ±-\/iV error bars), 
formed by binning the data Unearly in logE and 
normaUzing, is compared with the prediction (1) 
for lognormal field statistics. The symbols and 
dashed line are as in Figure 2. 



Fig. 4. — The distribution P{X) for all data in 
phases 485 - 500 (crosses), inclusive, is compared 
with the SGT prediction (solid line). See text for 
details. 
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